What should a statistical mechanics satisfy to reflect nature? * 
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There is no compelling reason imposing that the methods of statistical mechanics should be 
restricted to the dynamical systems which follow the usual Boltzmann-Gibbs prescriptions. More 
specifically, ubiquitous natural and artificial systems exhibit complex dynamics, for instance, generic 
stationary states which are not ergodic nor close to it, in any geometrically simple subset of the a 
priori allowed phase space, in any (even extended) trivial sense. A vast class of such systems appears, 
nevertheless, to be tractable within thermostatistical methods completely analogous to the usual 
ones. The question posed in the title arises then naturally. Some answer to this complex question 
is advanced in the present review of nonextensive statistical mechanics and its recent connections. 



A basic answer to the question of the title could of 
course be "to make predictions that are experimentally 
confirmed" . Still, what we would really like to know is 
what a priori mathematical properties or features should 
a formalism satisfy in order to enable useful statistical 
predictions in a manner similar to that of Boltzmann- 
Gibbs (BG) statistical mechanics. This task is either 
impossible or of extremely high difficulty. In order to 
be able to formulate an answer which could at least be 
partially satisfactory, let us analyze the basic structure 
of the BG theory. Indeed, no logical-deductive path ever 
existed for proposing a new physical theory, or for gen- 
eralizing a pre-existing one. In fact, such proposal fre- 
quently — perhaps always — occurs on a metaphorical 
basis. Therefore, the analysis of the structure of the BG 
theory will hopefully provide us a metaphor for formulat- 
ing a statistical mechanics which could perhaps be more 
powerful than the one we already have thanks to the ge- 
nius of Boltzmann, Gibbs, and others. It is perhaps wor- 
thy at this stage to explicitly declare that we are talking 
of a generalization of the BG theory, hy no means of an 
alternative to it. 



I. LEARNING FROM BOLTZMANN-GIBBS 
STATISTICAL MECHANICS 



A. Differential equation 

Which is the simplest ordinary differential equation? 
No doubt it is 



dx 



, 



(1) 



dx 



whose solution is y 
tempting to say 



^ = 1, (2) 



And the next one? It is 



dx ' 



(3) 



whose solution la y — e^ . Its inverse is y = Inx, which 
coincides, by the way, with the celebrated Boltzmann 
formula (called Boltzmann principle by Einstein) 



Sbg = khiW 



(4) 



where k is Boltzmann constant, and W is the measure of 
the space where the system is allowed to "live", taking 
into account total energy and similar constraints. So, 
typically, if we have an isolated TV-body Hamiltonian 
system (microcanonical ensemble in Gibbs notation), 
W is the dimensionless Euclidean measure (i.e., (hy- 
per)volume) of the fixed-energy Riemann (hyper) surface 
in phase space (Gibbs' F-space) if the system microscop- 
ically follows classical dynamics, and it is the dimension 
of the associated Hilbert space if the system microscop- 
ically follows quantum dynamics. In what follows we in- 
distinctively refer to classical or quantum systems. We 
shall nevertheless use, for simplicity, the wording "phase 
space" although we shall write down formulas where W 
is a natural number. 

If we introduce a natural scaling for x (i.e., if x carries 
physical dimensions) we must consider, instead of Eq. 
(3), 



whose solution (with ?/(0) = 1) is ?/ = 1. What could be 
considered as the second in simplicity? Surely 



dy 
dx 



(5) 
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in such a way that ax is a dimensionless variable. The 
solution is now 



and 



y 



(6) 



This differential equation and its solution appear to ad- 
mit at least three physical interpretations that are crucial 
in BG statistical mechanics. Let us state them without 
proof. The first one is {x, y, a) {t, ^, A), hence 



(7) 



where t is time, ^ = lim^x(o)^o ax{o) sensitivity 
to initial conditions, and A is the (maximal) Lyapunov 
exponent associated with a typical phase-space variable 
X (the dynamically most unstable one, in fact). This 
sensitivity to initial conditions (with A > 0) is of course 
the cause of the mixing in phase space which will guar- 
antee ergodicity, the well known dynamical justification 
for the entropy (4). 

The second physical interpretation is given by 
{x,y,a) — > (f, il, — 1/r), hence 



-t/T 



(8) 



where = 0(0] ; and r is the characteristic time 

associated with the relaxation of a typical macroscopic 
observable O towards its value at the possible stationary 
state (thermal equilibrium for BG statistical mechanics). 
This relaxation occurs precisely because of the sensitivity 
to initial conditions, which guarantees strong chaos (es- 
sentially Boltzmann's 1872 molecular chaos hypothesis). 
It was apparently Krylov the first to realize [1] , over half 
a century ago, this deep connection. Indeed, r typically 
scales like 1/A. 

The third physical interpretation is given by (a;, y, a) — > 
{Ei,Zpi,-P), hence 



-l3Ei 



(9) 



where Ei is the eigenvalue of the i-th quantum state 

of the Hamiltonian (with its associated boimdary condi- 
tions). Pi is the probability of occurrence of the i-th state 
when the system is at its macroscopic stationary state 
in equilibrium with a thermostat whose temperature is 
T = l/k(3 (canonical ensemble in Gibbs notation). It is 
a remarkable fact that the exponential functional form of 
the distribution which optimizes the BG generic entropy 



w 



Sbg = -k^pilnpi , 

i=l 

with the constraints 

w 



(10) 



(11) 



w 



^^PiEi = U {U = internal energy), (12) 

i=l 

precisely is the inverse functional form of the same en- 
tropy under the hypothesis of equal probabilities, i.e., 
Pi = l/iy(Vi), hence the logarithmic Eq. (4). To the best 
of our knowledge, there is (yet) no clear generic mathe- 
matical linking for this fact, but it is nevertheless true. 
It might seem at first glance a quite bizarre thing to do 
that of connecting the standard BG exponential weight 
to the solution of a (linear) differential equation, in con- 
trast with the familiar procedure consisting in extremiz- 
ing an entropic functional (Eq. (10)) under appropriate 
constraints (Eqs. (11) and (12)). It might be helpful 
to remind to those readers who so think that it is pre- 
cisely through a differential equation that Planck heuris- 
tically found the celebrated black-body radiation law in 
his October 1900 paper [2], considered by many as the 
beginning of the path that led to quantum mechanics. 

Let us conclude the present remarks by saying that, 
when we stress that Eqs. (7), (8) and (9) naturally co- 
emerge within BG statistical mechanics, we only refer to 
the generic (or more typical) situations, not to all the 
situations. It is known, for example, that relaxation oc- 
curs through a power-law function of time at any typical 
second-order phase transition, whereas the Boltzmann- 
Gibbs weight remains exponential. 

B. Mean value 

The BG entropy (10) can be rewritten as the following 
mean value: 



Sbg = A:(ln— ) , 
Pi 



(13) 



w 



where (• • •) = ^^ft(- • •)• The quantity ln(l/pj) is some 

i=l 

times called surprise [3] or unexpectedness [4]. We notice 
that the averaged quantity has the same functional form 
as that corresponding to the equal probability case Eq. 
(4), where 1/pi plays the role of W. 

C. Entropy composition law for independent 
systems 

Let us consider systems A and B as probabilistically 
independent, i.e., such that P^'^^ = pfpf (V(i,j)). We 
can immediately prove that entropy (10) satisfies the fol- 
lowing property 



SsciA + B) = SsciA) + SsGiB) 



(14) 



i=l 



referred from now on as extensivity. This property is 
sometimes referred to as additivity, reserving the word 
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extensivity for the infinitely many body systems; we will 
for simplicity not make such distinction here. 

The imear property (14) of course encompasses the fact 
that, since Wa+b = WaWb, whenever we have equal 
probabilities, the logarithmic form (4) is absolutely fit- 
ted. For example, if wc have N independent coins (or 
dices), it is = 2^ (or 6^), hence Sbg = Nk\n2 
(or Sbg = Nkln6). If we have, as another example, a 
d = 3 regular lattice with ferromagnetic Hciscnbcrg inter- 
actions between first neighbors at very high temperature, 
it is ~ Ap^ (with ^ > 0, p> 1, and N oo), hence 
Sbg ^ Nklnp. In all these cases, we have Sbg N, 
which precisely fits the Clausius concept of thermody- 
namic entropy. We shall discuss in Section II what can 
we do when the ubiquitous (but not obligatory) behavior 
W{N) ~ (with N » 1) is drastically violated, e.g., 
when W oc A'"'', with 7 > 0, which also appears to be 
ubiquitous in both natural and artificial systems. 

D. Concavity 

Let us consider two probability distributions {pi} and 
{p'^ for a given system (z = 1, 2, • • • , W). We shall define 

an intermediate distribution as follows: 



Pi = m + (1 - fJ')Pi (0 < At < 1) 



(15) 



An entropic functional S{{pi}) is said concave if and only 
if 

s{{Pi}) > + (1 - M)^({p-}) (v{ft},v{pa,VM) . 

(16) 

It can be easily shown that Sbg is concave. This property 
generically yields, within BG statistical mechanics, ther- 
modynamical stability, i.e., stability of the system with 
regard to energetic perturbations. For the canonical en- 
semble this implies S bg / dU^ < (i.e., positive specific 
heat). For the microcanonical ensemble, the situation 
can in fact be more complex (see, for instance, [5]). Also, 
it is this property which makes that if we put into ther- 
mal contact two systems which are initially at different 
temperatures (of the same sign, usually positive), ther- 
mal equilibrium occurs at a temperature which is nec- 
essarily intermediate between the two initial ones. And 
this occurs with a variation of the total entropy which is 
necessarily an increase. We therefore see how central the 
concavity of the entropy is for both the 0*'* and the 2"'^ 
principles of thermodynamics. 

E. Stability or experimental robustness 

Lcsche [6] addressed in 1982 a very basic (and, cu- 
riously enough, widely unknown) property. A necessary 
condition for a (positive) statistical functional 0({pj}) to 
be a physical quantity is [6] that, under arbitrary small 
variations of the probabilities {pi}, its relative variation 
remains small. We then say that 0{{pi}) is stable [6] (we 



shall hereafter also use the expression experimentally ro- 
bust, in order to avoid confusion with the thermodynamic 
stability mentioned in the previous subsection). To be 
more precise, if we consider two probability sets {pi} and 
{p-} associated with W microstates, the measure of the 
size of the deformation can be defined as follows [6] 



\\P-P'\\=^\pi-P'i 



(17) 



The condition of stability or experimental robustness of 
0{{pi}) is then given by 



\\P-P'\\<5e 



R: 



0{{Pi})-0{{p'^) 



<e , 
(18) 



for any £ > 0, with 6s > being independent from 
W, Omax being the maximal value that 0{{pi}) can at- 
tain. This implies, in particular, that lim lim R = 

E— >0 W^OO 

lim lim i? = 0. 

It can be shown [6] that Sbg (for which Omax = 
k\nW) is experimentally robust. 

F. Finite entropy production per unit time 

Let us consider the allowed phase space (assumed con- 
tinuous and finite) of our system, and let us partition 
it in >> 1 little cells. Assume next an ensemble of 
M >> 1 initial conditions, all inside one of the W little 
cells, and let the microscopic dynamics of the system to 
run. If the system is strongly chaotic (i.e., if it has posi- 
tive Lyapunov exponents) , the points will quickly spread 
everywhere inside the allowed region. We might define 
a set of probabilities pi = Mi{t)/M {i = 1,2, ■■■,]¥), 
where Mi{t) is the number of points inside the i-th cell 

at time t ^^Mi(t) = M, \ftj . The entropy SBG{t) will 

be zero at t = 0, and will start increasing as t goes on, 
finally saturating at some value (which is klnW if the 
system is ergodic; more precisely, if it equally visits all 
regions of the allowed phase space). The quantity 



Kbg = lim lim lim 

t— >cx) M— >o< 



SBG{{Pi{t)})/k 
t 



(19) 



will in general be finite and equal to the Kolmogorov- 
Sinai entropy (which constitutes a slightly different, more 
familiar in the mathematical community, definition of the 
same concept). This finitcncss follows from the Pcsin the- 
orem, which can be rewritten essentially as follows 



Kbg = ^2 ' 



(20) 



where r runs over all the positive Lyapunov exponents of 
the system. 
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II. USING WHAT WE LEARNT FROM 
BOLTZMANN-GIBBS STATISTICAL 
MECHANICS TO GENERALIZE IT 

There are several other properties than those discussed 
above, which also specifically characterize BG statistical 
mechanics, but we shall restrict the present analysis to 
only those, i.e., differential equations, mean value, en- 
tropy composition law, concavity, experimental robust- 
ness, and entropy production. For further — and con- 
sistent — characterizations, sec [7]. As we already men- 
tioned, there is of course no logical-deductive manner to 
generalize a physical theory. Or, to put it on more general 
grounds, there is no generic or unique way of generaliz- 
ing a logically consistent set of axioms into another one 
which also is logically consistent and which, by construc- 
tion, recovers the original one as a particular case. It is 
therefore only metaphorically that we shall use, in what 
follows, the mathematical structure of BG statistical me- 
chanics in order to generalize it. 

A. Differential equations 

Can we unify Eqs. (1), (2) and (3) in a single differen- 
tial equation? Yes, through 



dy 
dx 



= a + by . 



(21) 



Can we do it minimally, with only one parameter? (in- 
stead of two, namely a and b). Yes, out of linearity, 
through 



(22) 



Eqs. (1), (2) and (3) are respectively recovered for 
q — oc, g = and g = 1. The solution of Eq. (22) 
(with 2/(0) = 1) is given by 



y = [1 + (1 - q)x] 



1/(1-9) 



< (e? = e-). 



(23) 



The inverse function of the q-exponential is the q- 
logarithm, defined as follows 



y = 



!;l-9 _ 1 

1-a 



= lug x (Ini X = hix) . 



(24) 



It should be clear that these arguments about differen- 
tial equations have by no means a provatory nature nor 
intention. They are given here to provide some specific 
"feeling" about linearity and nonlinearity, thus providing 
some intuitive plausibility to the generalization that we 
propose in what follows. On these grounds, we may ex- 
pect the Boltzmann principle, Eq. (4), to be generalized, 
for equal probabilities, as follows 

Sgipi = 1/W, Vi) = klngW = fc ^j^~-^ , (25) 



As for the BG case, if x carries a physical dimension, 
we must consider, instead of Eq. (22), 



hence 



dy g 

Tx = 



y 



(ai 



(26) 



(27) 



As for the BG case, we expect this solution to admit at 
least three different physical interpretations. The first 
one corresponds to the sensitivity to initial conditions 



(28) 



where \q generalizes the Lyapunov exponent or coeffi- 
cient. Expression (28) was conjectured in 1997 [8], and, 
for unimodal maps, proved recently [9,10]. 

The second interpretation corresponds to relaxation, 



I.e., 



(29) 



There is (yet) no proof of this property, but there are 
several verifications (see, for instance, [11] for a quantum 
chaotic system). 

The third interpretation corresponds to the energy dis- 
tribution at the stationary state, i.e., 



-/9,-Ei 



Pi = 



W 

j=i 



(30) 



This is precisely the form that comes out from the op- 
timization of the generic entropy Sq under appropriate 
constraints [12-14]. This form has been observed in a 
large variety of situations (see [15] for mini-reviews). 

Before closing this subsection, let us stress that there 
is no reason for the values of q appearing in Eqs. (28), 
(29) and (30) be the same. Indeed, if we respectively de- 
note them by qsen {sen stands for sensitivity), qrei {rel 
stands for relaxation) and qstat {stat stands for station- 
ary state), we typically (but not necessarily) have that 
qsen < 1, qrei > 1 and qstat > 1- The possible con- 
nections between all these entropic indices are not (yet) 
known in general. However, for the edge of chaos of the 
z-logistic maps (see [16-19] and references therein) we 
do know some important properties. If we consider the 
multifractal f{a) function, the fractal or HausdorfF di- 
mension df corresponds to the maximal height of /(a); 
also, we may denote by amin and ctmax the values of a 
at which f{a) vanishes (with amin < ctmax), It has been 
proved [9,16] that 



1 



1 



1 - 



(31) 



Moreover, there is some numerical evidence [18] suggest- 
ing 
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— (X (1 - rf/) . 



(32) 



Unfortunately, we know not much about qstat, but it 
would not be surprising if it was closely related to qrel- 
They could even coincide, in fact. 

B. Mean value 

Since we have seen in the previous subsection that 
the logarithmic function naturally generalizes into the 
g-logarithmic one, let us define 



Pi 



(33) 



where we may call lng(l/p,) the q-surprise or q- 
unexpectedness. Then, using Eq.(24), it is straightfor- 
ward to obtain 



Sq — k 



w 

i=l 



1 



{Si = Sbg) 



(34) 



which is the entropy on which we shall base the present 
generalization of BG statistical mechanics [12-14] . This 
entropy of course recovers Eq. (25) for equal probabili- 
ties. 

C. Entropy composition law for independent 
systems 

If we consider now the same two probabilistically in- 
dependent systems A and B that we assumed before, we 
straightforwardly obtain 



SgjA + B) 

k 



(1-9) 



k 



(35) 



We re-obtain Eq. (14) in the limit {1 — q)/k ^ 0. Since 
Sq is always nonnegative, we see that, if q < 1 {q > 1), 
we have that Sq{A + B) > Sq{A) + Sq{B) {Sq{A + B) < 
Sq{A)+Sq{B)), which shall be referred as the superexten- 
sive (subextensive) case. It is from this property that the 
expression nonextensive statistical mechanics was coined. 

Let us address now the case W oc N'^{j > 0; >> 1) 
that we mentioned earlier. If we replace this into Eq. (25) 
we obtain Sq oc A^t(i~9) if g < 1. Consequently, it exists a 
special and unique value of q, namely q* = 1 — (I/7) < 1, 
for which we have Sq- oc N, once again in agreement 
with the Clausius idea of entropy! This only occurs be- 
cause the case W oc is incompatible with any hy- 
pothesis of probabilistic independence (even in the limit 
N ^ 00). Amusingly enough, we see that it might hap- 
pen that the entropy Sq with q ^ 1, which is nonexten- 
sive for independent systems, can be extensive when we 
consider a special class of systems which include strong 
dependence. More explicitly, if we asymptotically have 



W cx p^, it follows (for and Nb large and indepen- 
dent) W{A + B) oc p^^+^B = pNApNu ^ W{A)W{B), 
hence we must choose g* = 1 to have a finite limit 
for the entropy per particle. Whereas, if we asymp- 
totically have W oc A^^, it follows that W{A + B) oc 
{Na + N^Nl oc W{A)W{B), and the way for 

still having a finite limit for the entropy per particle is 
choosing q* = 1 — (I/7). 

It is interesting to notice [10] that if we replace W oc 
N'^ into S2-q we obtain S^-q oc N'^^'^~^\ hence S2-q** oc 
Nifq** =1 + (1/7) > 1. Naturally q** = 2- q* . 

In fact, this remark can be trivially generalized. Con- 
sider K(g), where k is any monotonic continuous function 
satisfying «;(1) = 1. If we focus on the decreasing ones, an 
example of such function is the one already considered, 
namely k = 2 — q. Another one is k = 1/g. Other possi- 
bilities can be obtained by successively applying the two 
just mentioned, i.e., k = 1/(2— 1/q) = q/{2q— 1), or, the 
other way around, k = 2 — 1/(2 — (?) = (3 — 2q)/{2 — q). 
(In fact, the four possibilities we have just considered 
belong to the class n = {a + hq)/[{a + 2h)q — 6] with 
(a, &) e V?, which is symmetric with regard to the axis 
K = q). Then, for the case W oc we are considering, 
we have that ^^(q*) oc N with q* ~ k^^{1 — (I/7)), 
denoting the inverse function. If K(q) is an increasing 
(decreasing) function, we have that q* < l{q* > 1). The 
considerations we have done in this paragraph might be 
not unrelated with the possible connections between the 
entropic indices qrei, qstat and qstat mentioned earlier. 

Finally, let us remind that Eq. (35) can be rewritten 
in the following extensive form [12]: 



S^{A + B) = S^{A) + S^{B) 



(36) 



where 



w 



ln[l + (1 - q)Sq/k] 
1-q 



1-q 



(37) 



is the Renyi entropy. This entropy is extremal &i 1/W 
(equal probabilities), and attains the value \nW for all 
values of q. Clearly, this fact makes it useless for having 
a finite entropy per particle in the case W oc . It is, 
nevertheless, an interesting functional for geometrically 
characterizing multifractals, as long known. 

D. Concavity 

It can be shown that Sq{{pi\) is a concave (convex) 
functional for all positive (negative) values of q. It also 
follows that, for the canonical ensemble, the specific heat 
is necessarily positive for g > 1, but not necessarily for 
< g < 1. 

Let us mention that Renyi entropy is concave for g < 1 , 
but has no definite concavity (or convexity) for 5 > 1. We 
consider this as a serious drawback for the use of Renyi 
entropy with q > 1 for thermostatistical and thermody- 
namical purposes. 
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E. Stability or experimental robustness 



(39) 



Abe has recently proved [20] a remarkable fact, namely 
that Sq is stable (experimentally robust) for all positive 
values of q. In contrast, Renyi entropy is unstable (exper- 
imentally fragile) for all values oi q ^ \ [6] . Once again, 
we consider this a further serious drawback for the use of 
the Renyi functional as a basis for thermodynamics. 

F. Finite entropy production per unit time 

If the dynamics of the system is such that the Lya- 
punov exponents vanish, Eq. (20) provides = 0, which 
clearly is poorly informative. We would like to "unfold" 
this trivial equality, and know more about those zeros. 
The functional Sq has been shown to enable precisely 
this in a variety of situations, including various low- and 
high-dimensional nonlinear dynamical systems, both con- 
servative and dissipative, both discrete and continuous in 
time and space. The situation can be summarized as fol- 
lows. 

We generalize definition (19) as follows: 



Ka 



hm hm hm M^lMlA 

t^oo W— >oo M^oo t 



K 



BG) 



(38) 



What has been verified for vast classes of Kbg = non- 
linear dynamical systems (see, for instance, [21]) is that 
a special and unique value of q exists (coincident with 
qsen < 1 every time checking has been performed) such 
that ifg,^„ is finite, whereas Kg vanishes (diverges) for 
all values of 5 > qsen (9 < qsen)- In other words, for such 
systems, the entropy whose production per unit time is 
finite is Sq^^^, not Sbg- 

The full characterization of the systems so behaving, 
and the full comprehension of their dynamical details in 
what concerns wandering in phase space, are still lack- 
ing. The scenario seems nevertheless relatively clear. If a 
chaotic sea (frequently single connected) exists in phase 
space, and we perform an average over many initial con- 
ditions (i.e., many choices for the little cell, among the W 
cells existing in the partition, where we initially place the 
M points) all over the entire allowed phase space, then 
the qsen = 1 entropy production is finite, independently 
on whether the m,easure of occupancy is uniform (i.e., er- 
godic in the entire allowed phase space) or not. But if 
the evolution in phase space is such that the system re- 
mains sensibly long times in regions where the structure 
is multifractal-like (possibly scale-free-like, in the sense 
of Barabasi et al, see [22] and references therein), then 
one expects qsen < 1- More details will be shown later in 
this paper. 

The generalization of the Pesin theorem (Eq. (20)) 
along the present lines was conjectured in 1997 [8], and 
has been recently proved by Baldovin and Robledo [10] 
for unimodal one-dimensional maps. More precisely, they 
have proved that 



We believe this to be one of the cornerstones of the entire 
theory. 

It is clear that the comments we did earlier concerning 

the behavior of Sq as a function of N are applicable here 
as a function of t. In other words, the entropy produc- 
tion per unit time is finite also for S'K(g) (t), i.e., if Kq^^^ is 
finite, the same occurs for K^^q^^^y Also, it is clear from 
Eq. (37) that for all those g 7^ 1 systems for which Sq{t) 
is asymptotically linear, Renyi entropy is not. There- 
fore, for such systems, there is no value of q (excepting of 
course whenever qsen = 1) for which S^ asymptotically 
could yield a finite entropy production per unit time. 



III. THE CANONICAL ENSEMBLE AND 
CONNECTION TO THERMODYNAMICS 

In the previous Sections we have addressed the ques- 
tion of which one, among the infinite possible general- 
izations of SbGi we want to postulate as the basis for 
generalizing BG statistical mechanics. We have decided 
it will be Sq. Now, if we have an entropic form, in princi- 
ple any entropic form, and are interested in cybernetics, 
control theory, information theory, and related matters, 
there are many things that we can do just with that. 
If we are, however, physicists, we might naturally think 
of doing statistical mechanics and thermodynamics. In 
other words, we will consider the energy (which repre- 
sents the physical support of the system) in addition to 
the entropy (which represents our information about that 
physical support). To do so, we shall first address what 
we consider to be the "Sancta Sanctorum" of statistical 
mechanics: conservative Hamiltonian systems, i.e., the 
systems Boltzmann and Gibbs themselves had in mind. 
More specifically, we shall address a large system in con- 
tact with an even (much) larger thermostat [canonical 
ensemble). Following along Gibbs' variational path, we 
shall extremize Sq with the norm constraint (11), and 
with a supplementary constraint related to the Hamilto- 
nian H of the system, namely [14] 



{n)q^ 



EL. I 



(40) 



with (...)i = (...); Pi = pfj Y^=iP'j is referred to as the 
escort distribution [23], and {Ei} is the set of eigenvalues 
associated with Ti and the corresponding boundary con- 
ditions. We discuss below the reasons which make desir- 
able the use of Pi instead of just pi for defining the energy 
constraint. The extremizing distribution corresponds to 
(meta) equilibrium and is straightforwardly shown to be 
given by 



-/3,(Ei-i7,) 



(41) 
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with 



3=1 



and 







(42) 



(43) 



P = 1/kT being the Lagrange parameter associated with 
constraint (40). We easily verify that q = 1 recovers the 
standard BG weight, g > 1 imphes in a power-law tail at 
high values of Ei, and Q' < 1 implies in a cutoff sA, high 
values of Ei. The (met a) equilibrium distribution (41) 
can be rewritten as follows: 



with 



w 



and 



0'a 



1 + (1 - q)pqUg 



(44) 



(45) 



(46) 



This form is particularly convenient for many applica- 
tions where comparison with experimental or computa- 
tional data is involved. 

Prom the preceding results the connection to ther- 
modynamics can be derived. In fact, the entire Lcg- 
endre transformation structure of thermodynamics is q- 
invariant. In particular, it can be proved that 



as well as 



The form adopted for constraint (40) (instead of the 
usual Eq. (12)) is at first sight astonishing, and surely 
demands clarification. It satisfies a remarkable set of (in- 
tertwined) properties, which we list now. 

(i) It satisfies, exactly as it happens when using pi, 
the basic property that the mean value of a constant is 
the same constant. This was not so in the version devel- 
oped in 1991 [13]. Chronologically speaking, to satisfy 
this property only became compelling when it became 
clear that the ^-generalization concerned a problem of 
lack of ergodicity, and not any unusual norm-preservation 
anomaly. 

(ii) It makes that the addition of microscopic ener- 
gies of independent systems preserves, at the macroscopic 
level, exactly the same form. In other words, Ef^'^^ = 

Ef + Ef with pf+^ = pfpf implies U^+^ = C// + . 
This property surely has an important role for the present 

generalization to satisfy the first principle of thermody- 
namics. This property was not satisfied in the 1991 ver- 
sion [13]. 

(iii) The distribution (41) obtained for the stationary 
state is invariant under change of the zero of energies, 
i.e., under uniform translation of the microscopic ener- 
gies. In more specific terms, if we add Eq to all energies 
Ei, the same Eq is, through Eq. (40), added to Uq, hence 
Ei — Uq docs not change, nor does pi. This property was 
not explicitly satisfied in the early 1988 and 1991 versions 
[12,13]. 

(iv) Abe and Rajagopal showed [24] that the tradi- 
tional steepest descent method (long ago used by Dar- 
win and Fowler to discuss BG statistics) naturally leads 
to the escort distribution in the energy constraint. It is 
basically related to the simple property de^/dx = (e^)'. 

(v) The form of constraint (40) makes that, under the 
optimization of Sq, the Lagrange a-parameter (the one 
associated with the norm constraint) appears in a func- 
tion which naturally factorizes out of the sum over all 
states {i}. This property warranties the definition of a 
partition function which, as usually, depends on the La- 
grange /3-parametcr but not on a. This interesting math- 
ematical fact was present in the 1991 version, but not in 
the version developed in 1988. 

(vi) If we express the two q > I canonical-ensemble 
constraints in the space of the energy, for high energy 
they asymptotically behave as indicated now: 



where 



ln„ Z„ = lUg Zg - pUq . 



Also, it can be proved that 



Y,Pi- dEg{E)p{E) 

J constant 



OC 



/ dEg{E)/E^/^'J-'^^ 

J constant 



as well as 



= ^dSq_ _ dUq 



QT2 • 



and 



J2Pi^i/J2P> dEg{E)E/E'i'^'i-^^ 

J constant 

* J 

= / dEg{E)/E' 

J constant 



^1/(9-1) 
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where g{E) is the density of states. Wc can observe 
that the domain of q where the two constraints are finite 
(hence mathematically and experimentally well defined) 
is the same; such nice property does not occur if we define 
the energy constraint using pi instead of the escort distri- 
bution. For example, if we assume the simple case where 
the high energy approximation is given by g{E) ~ 
{5 e SR), the convergence of both integrals is guaranteed 
for l/{q- 1) - S > 1, hence q < {2 + S)/{1 + 5). This 
feature constitutes an ingredient of consistency within 
the theory, which was not present in the 1988 version. 
Physically speaking, the internal energy Uq can be seen 
as a measure, at a given temperature, of the width or 
spread of the distribution pi above the lowest possible 
energy. If we show to a practically-minded scientist say 
an exponential distribution [p{x) oc e"""^ with a > for 
X > 0, and p{x) ~ for x < 0] and a power-law distri- 
bution oc 1/(1 + hx) with 6 > for x > 0, and 
p{x) = for X < 0], and ask him (her) what roughly 
arc the widths, he (she) will promptly check the width 
at about half value of the maximum for both cases, quite 
independently from the fact that the first moment of the 
first example is finite, whereas it diverges for the second 
one. The expectation value of the energy calculated with 
the escort distribution Pi, and not with pi, precisely is a 
measure for such width, for all cases, independently from 
what happens for the first moment of the distribution. 
In other words, this is a robust manner for characterizing 
this particular constraint. 

(vii) Let us analyze the same property as in (vi) but 
for a case where (x) =0, and consequently we need to 
refer to the second moment (x^). The already mentioned 
Levy-like superdiffusion illustrates such a situation. If we 
show to the same practically-minded scientist as before 
say a Gaussian [p(x) oc e~°^ ] and a Lorentzian [p(x) oc 
1/(1 -|- 6x^)] distribution, and once again ask him(her) 
what roughly arc the widths, he (she) will once again 
check the width at about half value of the maximum for 
both cases, quite independently from the fact that the sec- 
ond moment of the first one is finite, whereas it diverges 
for the second one. Fixing, for the present Levy-like ex- 



ample. 



precisely 



jr^dx x^pqix)]^\/[jr^dx [pq{x)r_ 

characterizes the width for all values of q below 3 (which 
simultaneously is the upper bound for normalizability!) 
[25-27] . The Lorentzian example just evoked corresponds 

to g = 2. 

All the above (i-vii) points were satisfactorily settled 
in the 1998 version of nonextensive statistical mechan- 
ics [14], which is the one that wc use since then, hence 
in the present review. For clarity on the situation, see 
the Table. However, in spite of the remarkable mathe- 
matical consistency of the 1998 version that one might 
check in this Table, the ultimate geometrical-dynamical 
understanding of escort distributions is still lacking and 
is object of current studies (see, for instance, [28]). The 
situation seems to have some analogy with the follow- 
ing gedanken problem. Let us consider a 10 cm linear 



object, and let us implement on it the construction of 
the well known Cantor set, where, at each hierarchy, 
the central third is deleted. The final object is a frac- 
tal with zero one-dimensional Lebesgue measure, and 
whose Hausdorff dimension is df = In 2/ In 3. What 
is the measure of such an object? It is clearly given 
by (10 cmYf = 10in2/in3 ^^in2/in3 ^ 427 cm°e^. 

In other words, we "start" with the entire Lebesgue 
measure 10 cm, but then we "correct" ! The function 
lUg W = {W^~'^ — 1)/(1 — g) does something quite analo- 
gous: it ""starts with the entire measure W of the Gibbs 
r space, and then it "corrects" . This would correspond to 
the fact that, for nonextensive systems, the entire phase 
space is dynamically not nearly entirely occupied (i.e., 
the system is not ergodic), but only a scale- free- like part 
of it (which depends on the initial conditions) is visited. 
Consistently, one expects some type of "correction". It 
is not impossible that the escort distribution would pre- 
cisely provide the necessary "correction". It is perhaps 
not unworthy to warn the reader that this is but a pos- 
sible scenario, which remains to be proved. 



IV. CALCULATING THE INDEX q A PRIORI 

For the present theory to be complete, it is obviously 
necessary to be able to calculate q a priori, i.e., from first 
principles. Consistently with the thought of Einstein [29], 
Cohen [30], Baranger [31], and many others, such a ba- 
sic and crucial task can be accomplished nowhere else 
than in the analysis of the microscopic dynamics (clas- 
sical, quantum, or others) of the system. However, very 
helpful connections can also be established at the level of 
Langevin-like, Fokker-Planck-likc, and similar equations, 
i.e., at the level of mesoscopic dynamics. In many oc- 
casions, natural and artificial systems whose microscopic 
or mesoscopic dynamics are either unknown or extremely 
complex have been addressed. For such systems, q has 
been obtained through the fitting of experimental data. 
We are not addressing these cases here, but more details 
on several such examples can be found in [15]. Our pur- 
pose here is to briefly mention several microscopic and 
mesoscopic determinations of q that are available in the 
literature. 

A. Determination of q from microscopic 
dynemiics 

(i) Low- dimensional dissipative systems (one- and two- 
dimensional dissipative maps) 

Let us consider the following one-dimensional dissipa- 
tive maps ( [16,17,32] and references therein): 

xt+i = 1 - a\xt\'' {z>l) , (47) 

referred to as the z-logistic map, and 

xt+i = dcos(7r]xt - l/2]^/2) (^ > i) . (48) 
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referred to as the z -periodic map. They belong to the 
same Feigenbauni-CouUet-Tresser universahty classes. 
More specifically, they share at their edge of chaos (e.g., 
for z = 2, Uc = 1.4011..., and dc = 0.8655...), the same 
Feigenbaum universal constants 6f{z) and apiz). 

Let us also consider ( [32] and references therein) a 
family of one-dimensional dissipative maps which belong 
to universality classes that are different from the just 
mentioned: 

Ot+i = O + [6»t - (l/27r)sm(27r(9t)]^/^ > 0; < Q < 1) . 

(49) 

referred to as the z -circle map. The edge of chaos oc- 
curs, for z = 3 (the usual case), at ilc = 0.6066.... For 
all three maps (47), (48) and (49), the following relations 
have been verified 

1 _ 1 _ 1 _ {z-l)hiaF{z) 

(z) In 6 

(50) 

where, as mentioned earlier, amin {otmax) is the minimal 
(maximal) value of a in the /(a) multifractal function 
(./(ttmin) = fioimax) = 0), and whcrc 6 = 2 for maps 
(47) and (48) and b = {^5 + l)/2 = 1.6180... {golden 
mean) for maps (49). 
It is clear that the dissipative Henon map 

xt+i = I - a\xt\'' + yt (51) 
yt+i = bxt (0 < \b\ < 1) 

belongs to the same universality class as the logistic map. 
Checking this case is not really necessary; nevertheless, 
the numerics have also been directly performed [33] . Con- 
sistency with the logistic map has indeed been verified. 

To be more explicit, qsen{z) has been calculated 
through various independent procedures, namely the sen- 
sitivity to the initial conditions [8,9], the multifractal 
function f{a) [16], and the finiteness of the entropy pro- 
duction per unit time [10,21]. A fourth method has 
been used as well, which we describe now. By measur- 
ing the shrinking of the Lebesgue measure, Moura et al 
obtained qrei{z) for the z-logistic maps. For example 
qrei{2) — 2.41. A connection has been recently obtained 
[19] between qsen and qrei, namely 

qreiiz, oo) - qreliz, W) oc 1/W^l9«"(-)I (52) 

where W is the number of cells in which the phase space 
has been partitioned (the larger W, the thinner the grain- 
ing). This equation can be used to simultaneously deter- 
mine qsen{z) and qrei{z) = qrei{z, oo). 

To illustrate all these methods, let us mention that, 
since api^) is known with at least 1018 digits for the 
logistic map, the corresponding value for 

9«en(2) = 0.2445... 



is known with the same number of digits. Cases like this 
one completely disqualify the critical — and, as we see, 
unfounded — remark that we hear occasionally, namely 
that q is no more than a "fitting parameter" . 

(a) Low- dimensional conservative systems (two- and 
four- dimensional conservative maps) 

The standard map is defined as follows: 

9{t-\-l) =p{t) + ^sm[2'Ke{t)]-\-eit) (modi), (53) 

p{t +1)= p{t) + ^ sin [2n9{t)] (mod 1) 

27r 

(a <E M, t = 0, 1, ...); 2ttp may in fact be regarded as 
the angular momentum of a free rotor subject to angle- 
dependent impulses of strength a at unit intervals of time. 

This two-dimensional map is conservative and simplec- 
tic. Also, it is symmetric with regard to p = 1/2. For 
|a| >> 1, the system is very chaotic (large Lyapunov ex- 
ponent) and ergodic. When a approaches Oc = 0.9716... 
from above, the region around p = 1/2 gradually be- 
comes isolated from the rest of the phase space. There- 
fore, if we define a kind of "dynamical temperature" 
T{t) = ((p(t))^) — (p{t))^, a nonzero-measure region ex- 
ists in phase space where, if we start inside it at i = 0, 
we observe two plateaux in T{t) [34]. The duration of 
the first plateau diverges when a ac. In other words, 
limt^oo lima^a^ r(a, t) \ima^aJ'^m^ocT{a,t). The 
existence of such a nonuniform convergence is, as we shall 
see later on, a distinctive feature which already suggests 
the possible existence of some degree of nonextensivity. 
Indeed, anomalies are observed in both the sensitivity to 
the initial conditions and in the entropy production per 
unit time. These anomalies appear to be consistent with 
each other, and they enable the characterization of an 
index qsen < 1- However, this index is only an effective 
one, and depends on a. For example, for a = 0.6, we can 
verify that qsen — 0.3, but a permanent drift of qsen is ob- 
served towards zero, when we average over many initial 
conditions all over the entire phase space and consider 
values of a gradually approaching zero. The limit a = 
corresponds of course to the integrable case of the map. 
Details can be found in [35]. 

The two-dimensional map has no Arnold diffusion. To 
verify the influence of dimensionality we also considered 
two standard maps coupled as follows: 

pi{t + 1) = pi{t) + ^ sin [2Trei {t)] , (54) 

ZTT 

(29 

P2{t + 1) ^ P2{t) + sin [27T02it)], 
ZTT 

diit + 1) = 9i{t)+pi{f + i) + h P2{t + 1), 

e2{t + 1) = 02(i) +P2{t + 1) + 6pi(i + 1), 

where ai,a2,b S M, t = 0,1,..., and all variables 
are defined mod 1. If the coupling constant b vanishes 
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the two standard maps decouple; if 5 = 2 the points 
i0i,Pi,e2,P2) = (0,1/2,0,1/2) and (1/2,1/2,1/2,1/2) 
are a 2-cycle for all (ai, 02), hence we preserve, in phase 
space, the same referential that we had for a single 
standard map. For a generic value of b, all relevant 
present results remain qualitatively the same. Also, we 
set ai = a2 = a so that the system is invariant under 
permutation 1 2. This four-dimensional map is con- 
servative and simploctic, and has Arnold diffusion. It 
exhibits two plateaus for T{t), like the two-dimensional 
case (with d^. = 0). The dependence of ?sen(^) is not 
very different from the two-dimensional case. We are 
presently studying iV >> 1 coupled maps to verify the 
effects of the thermodynamic limit. This might or might 
not stabilize a nontrivial value for qsen- 

Finally, let us mention in this section a connection 
with quantum chaos, namely for the quantum kicked 
top Hamiltonian. If we study the time evolution of the 
fidelity function (scalar product of the wave function 
of the original Hamiltonian with the wave function of a 
slightly perturbed Hamiltonian), three types of behaviors 
are typically observed, (i) In the regions corresponding 
to regular motion, the fidelity remains roughly constant; 
(ii) in the regions corresponding to standard chaos (with 
positive Lyapunov exponents for the associated classical 
system), the fidelity decreases exponentially with time 
(before entering the quantum interference region); (iii) 
at the border of the two regions, a power-law behavior 
has been recently observed [11] which precisely fits a q- 
exponential function. The value of qrei (associated with 
the zero perturbation limit) depends on the total mo- 
mentum value J of the top. It decreases with increasing 
J, and, in the limit J — > 00 (classical limit), qrei — 2.4 
can be extrapolated [11]. Further analysis is needed and 
welcome. 

(iii) High- dimensional dissipative systems (many-body 
dissipative models) 

A specific lattice Lotka-Volterra three-component 
model has been recently analyzed [36] in terms of entropy 
production per unit time. For d-dimensional growth, the 
following result has been obtained 

qsen = l-l/d. (55) 

Details for the g? = 1,2 cases can be found in [36]; the 
d = 3,4: cases have been studied by Anteneodo [37]. In 
fact the microscopic rules of the model are such that the 
growing droplet has a linear size which grows linearly 
with time t. The volume grows then as t''. If we apply to 
this quantity the 5-logarithmic function, we obtain that 
Sq scales like t if (1 — q)d = 1. From this observation, the 
result qsen = 1 — ^/d follows immediately [37] . The pos- 
sible connection with fractality as d varies is yet unclear 
[37]. 

Boltzmann d-dimensional Bravais lattice models for 
the Navier-Stokes equations for incompressible fluids 



have been recently revisited [38]. The imposition of the 
(physically desirable) Galilean invariance of the equa- 
tions mandates an unique entropy, and that entropy pre- 
cisely is Sq, with q depending on some details of the 
model {q is in fact analytically determined by a tran- 
scendental equation). For example, for the single- mass 
single-speed model, it is 

q=l-2/d. (56) 

(iv) High- dimensional conservative systems (many- 
body Hamiltonian systems) 

Hamiltonian systems is a central topic of statistical 
mechanics. We shall briefly review here some important 
results currently available for long-range-interacting clas- 
sical systems. Let us first focus on the inertia! XY fer- 
romagnetic model, characterized by the following Hamil- 
tonian [39,40]: 

i=i j 

where 9i is the i — th angle and pi the conjugate variable 
representing the angular momentum (or the rotational 
velocity since, without loss of generality, unit moment 
of inertia is assumed). Notice that the summation in 
the potential is extended to all couples of spins (counted 
only once) and not restricted to first neighbors; for d = 1, 
rij = 1,2,3,...; for d = 2, rij ~ l,\/2, 2,...; for d = 3, 
Vij = 1,^2,^3,2,.... The first-neighbor coupling con- 
stant has been assumed, without loss of generality, to be 
equal to unity. This model is an inertial version of the 
well known XY ferromagnet. Although it does not make 
any relevant difference, we shall assume periodic bound- 
ary conditions, the distance to be considered between a 
given pair of sites being the smallest one through the 2d 
possibilities introduced by the periodicity of the lattice. 
Notice that the two-body potential term has been written 
in such a way as to have zero energy for the global funda- 
mental state (corresponding to = 0, V?, and all 9i equal 
among them, and equal to say zero). The a ^ 00 limit 
corresponds to only first-neighbor interactions, whereas 
the a = limit corresponds to infinite-range interactions 
(a typical Mean Field situation, frequently referred to as 
the HMF model [39]). 

In the limit N ^ 00, the quantity TV = X^j^j 
converges to a finite value if a/d > 1, and diverges like 
j^i-a/d if < a/rf < 1 (like InN for a/d = 1). In other 
words, the energy is extensive for a/d > 1 and nonex- 
tensive otherwise. In the extensive case (here referred to 
as short range interactions; also referred to as integrahle 
interactions in the literature), the thermal equilibrium 
(stationary state attained in the t ^ 00 limit) is known 
to be the BG one (sec [41]). The situation is much more 
subtle in the nonextensive case {long range interactions). 
It is this situation that we focus on here. In order to 
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conform to the most usual writing, we shall from now on 
replace the Hamiltonian H by the following rescaled one: 



N 



n' 



(a > 0), (58) 



The molecular dynamical results associated with this 
Hamiltonian (now artificially transformed into an exten- 
sive one for all values oi a/d) can be trivially transformed 
into those associated with Hamiltonian H by re-scaling 
time (see [40]). 

Hamiltonian (58) exhibits in the microcanonical case 
(isolated system at fixed total energy U) a second order 
phase transition a,t u = U/N = 0.75. It has anomalies 
both above and below this critical point. 

Above the critical point it has a Lyapunov spectrum 
which, in the N ^ oo limit, approaches, for < a/d < 1, 
zero as N~'^, where K{a/d) decreases from 1/3 to zero 
when a/d increases from zero to unity, and remains zero 
for a/d > 1 [40,42]. It has a Maxwcllian distribution of 
velocities [43], and exhibits no aging [44]. Although it 
has no aging, the typical correlation functions depend on 
time as a g-exponential. Diffusion is shown to be of the 
normal type. 

Below the critical point (e.g., u = 0.69), for a nonzero- 
measure class of initial conditions, a longstanding qua- 

sistationary (or mctastable) state precedes the arrival 
to the BG thermal equilibrium state. The duration of 
this quasistationary state appears to diverge with N like 
N [43,45]. During this anomalous state, there is aging 
(the correlation functions being well reproduced by q- 
exponentials once again), and the velocity distribution is 
not Maxwellian, but rather approaches a g-exponential 
function (with a cutoff at high velocities, as expected for 
any microcanonical system). Anomalous superdifFusion 
is shown to exist in this state. The mean kinetic energy 
(oc T, where T is referred to as the dynamical temper- 
ature) slowly approaches the BG value from below, the 
relaxation function being once again a g-exponential one. 
During the anomalous aging state, the zeroth principle 
of thermodynamics and the basic laws of thermometry 
have been shown to hold as usual [46,47]. The fact that 
such basic principles are preserved constitutes a major 
feature, pointing towards the applicability of thcrmosta- 
tistical arguments and methods to this highly nontrivial 
quasistationary state. 

Although none of the above indications constitutes a 
proof that this long-range system obeys, in one way or 
another, nonextcnsive statistical mechanics, the set of 
so many consistent evidences may be considered as a 
very strong suggestion that so it is. Anyhow, work is 
in progress to verify closely this tempting possibility. 

Similar observations are in progress for the Heisenberg 
version of the above Hamiltonian [48], as well as for a 
XY model including a local term which breaks the an- 
gular isotropy in such a way as to make the model to 
approach the Ising model [49] . 



Lennard- Jones small clusters (with N up to 14) have 
been numerically studied recently [50] . The distributions 
of the number of local minima of the potential energy 
with k neighboring saddle-points in the configurational 
phase space can, although not mentioned in the origi- 
nal paper [50], be quite well fitted with g-exponentials 
with q = 2. No explanation is still available for this sug- 
gestive fact. Qualitatively speaking, however, the fact 
that we are talking of very small clusters makes that, 
despite the fact that the Lennard-Jones interaction is 
not a long-range one thermodynamically speaking (since 
a/d = 6/3 > 1), all the atoms sensibly sec each other, 
therefore fulfilling a nonextcnsive scenario. 

Finally, as a last example of Hamiltonian systems 
where nonextcnsive features have been analytically ob- 
tained (or numerically observed) we would like to men- 
tion anomalous transport in an optical lattice. The dis- 
tribution of velocities in such system has been recently 
shown [51] to be, for a specific microscopic model, a q- 
distribution with 



1 + 



44Ef 



Uo 



(59) 



where Er and Uq are microscopic parameters of the 
quantum optical problem. 

(v) Many-body growth models 

A growth model including preferential attachment has 
been recently introduced by Albert and Barabasi [22] as 
a prototype of emergence of the ubiquitous scale- free net- 
works. At each time step, m new links are added with 
probability p, or m existing links are rewired with prob- 
ability r, or a new node with m links is added with prob- 
ability 1 — p — r; all linkings are done with probability 



n(fc, 



(ki + 1)/X],(^j + l)i where ki is the number 



of links of the i — th node. The exact stationary state 
distribution of the number k of links at each site can be 
written (although apparently not identified by the au- 



thors) as p{k) oc e, 



with 



q = 



2m(2 -r) + l-p- 



(60) 



ro(3 — 2r) + 1 — p — r ' 
fco > being an explicit function of {p,r,m). 

(vi) Diffusion in the diluted hypercube 

The diluted iV-dimensional hypercube has been con- 
sidered since long. Its percolation '"threshold" is given 
[52] hy pc = a + f cr^ + i5^3 _^ ^ _ ^^leie a = 1/{N - 1). 
Strictly speaking, the quantity Pc is an effective perco- 
lation threshold; indeed the exact percolation threshold 
clearly is pc = 0. Lemke and Almeida [53] have recently 
studied the model where a random walker diffuses in this 
structure. More specifically they studied the entropy pro- 
duction for N increasing up to 23, and a finite value was 
found only for g = 1 at p = 1 (full hypercube), and for 
q = 0.056 at p = Pc- 
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B. Determination of q from mesoscopic dynamics 

(i) Nonlinear Fokker-Planck equations (Correlated 
anomalous diffusion) 

In a variety of physical situations [26], it is appropri- 
ate to consider the foUowing nonhnear Fokker-Planck-hke 
equation (sometimes referred to as the Porous Medium 
Equation): 

dp(x,t) d , , ^d'^\p(x,t)]^ , ^, 

(61) 

If we assume that at t = we have the paradigmatic (and 
quite usual) distribution p{x,0) — d{x), it can be shown 
[26] that, for F{x) = fci - k2X (ki G TZ; k2 > 0) and all 
{x,t), the (stable) solution is given by 

p(a;,t)ae-''W[^-^^Wl' (9 < 3) , (62) 

where /3(t) and XM{t) are smooth explicit functions of t, 
and 



q = 2-u (g < 3) . 



(63) 



(a) Fractional- derivative Fokker-Planck equations 
(Levy anomalous diffusion) 

It is known [54] that p{x) oc e~^^^ {q < 3) optimizes 

^ l-Jdx\p{x)]'' 



1 



(64) 



under appropriate constraints. If we convolute N times 

p{x) {N 00), we approach a Gaussian distribution if 
q < 5/3 and a Levy L^^(x) one if 5/3 < g < 3. Levy 
distributions are the solutions of the following equation 
djA^ = i^^^jt^ (0 < 7L < 2). The index 7^ of 
these Levy distributions is related to q as follows: 



5=^^ (5/3<g<3). 

7f, + 1 



(65) 



In general, the Levy distribution and the g-Gaussian 
coincide only asymptotically (the tails). The unique 
exception is for q = 27^, = 2, in which case they are 

identical. 

(Hi) Nonlinear fractional- derivative Fokker-Planck 
equations 



The equation 

dp{x,t) ^^d^[p{x,t)r 



dt 



dx^ 



(u G TZ; J € TZ) . (66) 



is an interesting one. The 7 = 2 particular instance 
reproduces the situation addressed with Eq. (61); the 



V = 1 particular instance reproduces the partial differ- 
ential equation above Eq. (65). We do not know the 
solution for the generic case. However, we do know it for 
the (new) case y = (2 — 7)/(l -|- 7) [55]. It is given by 



Q = 



7 -h 3 _ 5 -h 2!/ 
7+1 ~ 3 



(67) 



with < 7, < 2. As for the Levy case, it is only 
asymptotically that the solution coincides with the q- 
Gaussian function. More, and more complex, situations 
along these lines can be found in [56]. 

(iv) Anomalous Langevin equations 

We may consider the following Langevin-like equation 
[27] 

X = --fx\x\'^'^'-'^^ + x\x\'-^^{t) -\- r]{t) {s > 0), (68) 

where ^{t) and r]{t) are independent and Gaussian- 
distributed zero-mean white noises, satisfying 



and 



(^(*)^(*')) = 2M5(t - t') (M > 0) (69) 



{7]{t)r]{t')) = 2A6{t - t') (A > 0) . (70) 



M and A stand for m,ultiplicative and additive respec- 
tively. If 7 > M(l — s), the distribution corresponding 
to the stationary state is given by 



P{x) 



(t/3) + M 1^123 



with 



Q = 



(7/s) + 3M 
(7/s)-|-M 



(71) 



(72) 



Notice that, interestingly enough, the case s = M = 1 
reproduces the results indicated in Eqs. (65) and (67). 

Another interesting Langevin-like case is that of di- 
chotomic colored noise addressed in [57]. One may con- 
sider the following stochastic differential equation: 



(73) 



where ^{t) is a dichotomous noise of values ±a, satisfying 
{^{t + t)0 = a^e-^^l^l . (74) 

The stationary state probability is, although not writ- 
ten in this way in the original paper [57], given by 
p{V) oc e-^^' with 



(75) 
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/? being given by some smooth function of the model 
parameters. 

(v) Phenomenological approach of pressure fluctua- 
tions in multiphase flow 

A granular system has been experimentally studied [58] 
which consists of bubbling fluidized beds set in motion 
by a vertical gas stream. In addition to the experiments 
themselves, a phenomenological approach has been ad- 
vanced. It provides, for the distribution of the pressure 
fluctuations, a g-Gaussian with 

1='^ + yZJJ^ (t>3/2), (76) 

where r is a phenomenological parameter of the bub- 
ble size distribution. We must warn the reader that the 
authors argue [58] that, from their standpoint, no deep 
connection exists with nonextensive statistical mechan- 
ics. We do not necessarily agree with their interpreta- 
tion. But it is not the aim of this short note to address 
this particular point in depth. We just refer the existence 
of this mathematical result in connection with those ex- 
periments. 

(vi) Scalar granular gases 

Baldassarri et al [59] have numerically and analytically 
discussed some isolated one-dimensional granular gases 
in a recent work. A restitution coefficient r e [0, 1) has 
been asstuned for two-body collisions. The velocity dis- 
tribution in the asymptotic cooling regime of a specific 
pseudo-Maxwell (Ulam's) model, has been shown to be 
a distribution, which, although the authors apparently 
have not identifled it as such, precisely is a g-Gaussian 
with q = 3/2. 



V. TOWARDS A GENERALIZED 
NONEXTENSIVE STATISTICAL MECHANICS 

A statistical mechanics can be seen as basically con- 
sisting in the appropriate choice of a physical entropic 
functional S{{pi}). Once this is done, we must consis- 
tently choose the constraints (X]i=i-Pi — 1? ^iid possibly 
others) to be used for describing various physical situa- 
tions (isolated, or in equilibrium with an infinitely large 
reservoir, or others). We must then extremize (typically 
maximize) the entropy under these constraints, in order 
to obtain the probability distribution corresponding to 
the stationary state, whenever it exists. Relaxation phe- 
nomena towards the distribution associated with the sta- 
tionary state(s) can be characterized through S{{pi{t)}), 
among others. The stationary state entropy is of course 
given by S {{pi{oo)}) . The experimental information can 
provide (direct or indirect) evidence on {pi{t)}, and on 
S{t) = S{pi{t)}). But establishing the functional S'({pj}) 



itself demands the knowledge of S{t) for almost all possi- 
ble trajectories. So, in practice, the functional S{{pi}) is 
not deduced or measured, but found heuristically. This 
functional is in fact a mathematical construct which, for 
a variety of macroscopic physical quantities, can essen- 
tially replace (as the genius of Boltzmann showed to us) 
the detailed knowledge implied at the microscopic dy- 
namical level. 

The microscopic dynamics of an isolated system might 
be such that, for practically all possible initial conditions 
and as time increases, it quickly evolves towards an uni- 
form occupation of phase space (either of the entire phase 
space, as for the N ^ oo d = 3 short-range Heisenberg 
ferromagnet above its critical total internal energy, or of 
a simple and well defined part of it, as for the same ferro- 
magnet below its critical internal energy, once a specific 
symmetry has been broken). The system is then said er- 
godic, and Boltzmann's basic equal-probability hypothe- 
sis indeed applies. The entropy to be used in such a case 
is widely known to be Sbg- 

What happens, however, if the system has more com- 
plex dynamics, such that it tends, when isolated, to oc- 
cupy phase space in a nontrivial nonuniform manner? 
This is where generalizations of BG statistical mechanics 
possibly become a natural way out. This imposes, in our 
opinion, the need of physical entropies different from the 
BG one. 

The nonextensive entropy Sq{{pi}) 

(i) satisfies (Vg > 0), among others, three mathemat- 
ical properties (namely concavity, experimental robust- 
ness, and finiteness of the entropy production per unit 
time) that are by no means trivial to satisfy (e.g., Renyi 
entropy satisfies none of these Vg > 0), 

(ii) yields a stationary distribution (the (/-exponential 
function, which is asymptotically a power-law) which is 
ubiquitously found in natural and artificial systems, and 

(iii) is in principle completely determined once its en- 
tropic index q is calculated from microscopic or meso- 
scopic dynamics. 

It is therefore allowed to consider Sq as a quite strong 
candidate for physically extending the applicability of the 
methods of BG statistical mechanics to a specific class of 
dynamical systems among those who do not comply with 
the usual ergodicity requirements. The rigorous charac- 
terization of this class is yet to be better understood, 
but we already dispose of a good basis for conjecturing 
that it concerns systems that mix slowly in phase space, 
yielding a long-standing (multi) fractal-like occupation of 
it (possibly with a scale-free structure, as it is nowadays 
called in the field of networks) . Such is the case of various 
dissipative low-dimensional nonlinear dynamical systems 
at their edge of chaos, conservative low-dimensional non 
linear dynamical systems close to the frontier between 
chaos and integrability, (possibly) long-range-interacting 
Hamiltonian systems, and others. 

Are we compelled to stop here, or can we hope to 
similarly cover even more complex dynamics, eventually 
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yielding stationary states (or sequences of quasistation- 
ary states) whose phase space structure is even more com- 
plex than the one just mentioned? The bottom line is 
that we see no basic reason for stopping. However, a 
highly nontrivial problem must be solved before further 
progressing along this line, namely the proposal of an 
appropriate entropy functional. There are at least two 
lines of thought (hereafter referred to as the differential 
equation path, and the superstatistics path) that converge 
onto a solution of such a task. Let us next briefly mention 
both. 

A. The differential-equation path 

We have seen that BG and nonextensive statistical me- 
chanics are deeply related to the differential equations 
(5) and (26) respectively. Can we unify them in such a 
way as to have a crossover between one and the other, 
depending on the value of the independent variable x? 
Yes, we can, and this has been advanced in 1999 [60] 
for re-association in folded proteins, and recently used 
for a phenomenological theory for the flux of cosmic rays 
[61]. The associated differential equation we propose is 
as follows: 



dy , 



ar)y'i (y(0) = l) 



hence 



y 



1 - 



Oi 



(77) 



(78) 



Eq. (77) recovers Eq. (26) if we consider ax = 0. It re- 
covers Eq. (5) if we consider = ai, and also if we just 
consider q = 1. The existence of this last (convenient) 
possibility is the rc;ason for which we presented the gen- 
eralization in the form of Eq. (77), rather than in the 
totally equivalent form || = aiy + a'^y''. An interesting 
situation occurs for the case q < 1 with < ai << a^, 
as well as for the case q > 1 with Of, << ai < 0. 
For < a; << Xcrossover = [1/[(1 - q)aq] we have 
(nonextensive statistical mechanics), whereas 



we have y ~ e^^^ (BG statistical 



y ( 

for X ^ ^ Xf^j^ossover 

mechanics). 

We can further generalize Eq. (77) in a natural way, 
namely through 



+ (0,-0,02/" (y(o) = i) 



(79) 



whose solution is a combination of hypergeometric func- 
tions (sec [60]). As before, an interesting situation oc- 
curs for the case q < q' < 1 with < a,' << a,, as 
well as for the case q > q' with a, << a,/ < 0. For 

l-g' l-q 

Q < X « Xcrossover = [(1 " Q')aq]^^/[(1 - gOvl'^ 

e„'' (g- nonextensive statistical mcchan- 
crossov&r we have y ^ eq'' {q'- 



we have y 

ics), whereas for x >> x, 
nonextensive statistical mechanics). 



Before closing this subsection, let us make a remark 
concerning the history of physics. It seems quite reason- 
able to discuss the time dependence of quantities like 
the sensitivity to initial conditions and the relaxation 
of physical quantities in terms of a differential equation. 
But to do so for the energy distribution associated with 
a stationary state seems, as earlier mentioned, rather 
amazing. To make such procedure somewhat more ac- 
ceptable it might be interesting at this point to review 
a few equations of Planck's October 1900 paper [2]. He 
writes the following two equations: 



U{P + U) ' 



and 



dS_ 
dU 



1 
T 



(80) 



(81) 



where S is the entropy, U the internal energy, T the ab- 
solute temperature, a and (3 two coefficients to be fixed. 
The first of these two equations was heuristic. If we re- 
place the second equation into the first one, we obtain 



hence 



_d_ 

dU 



dU 



[/(/? + U) 



d{l/T) a a 



(82) 



(83) 



It happens that this equation precisely is the {q',q) = 
(1,2) particular case of our Eq. (77). Planck, in fact, 
does not write our present Eqs. (82) and (83). After 
writing our present Eqs. (80) and (81), he jumps (taking 
into account the density of states) to the writing of the 
celebrated law for the radiation of the black body. Al- 
though not always noticed, the bosonic nature of photons 
(i.e., the special constraint implied by the symmetriza- 
tion of the quantum wavefunctions) makes that Planck's 
law does undergo through a crossover from a power law 
{q = 2) at low energies to an exponential law (g = 1) at 
high energies. 

B. The superstatistics path 

It was noticed in 2000 by Wilk and Wlodarczyk [62] 
that the g-exponential distribution can be written as a 
Laplace transform involving the BG exponential weight. 
More precisely, 



dl3e-^^*f,{p) 



(84) 



where fq{(3) is the (also called the Gamma) distribu- 
tion. They further noticed that 



(85) 
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where (• • •) = dj3 {• • •)fq{(3). This observation was 
further developed by Beck [63] . The idea that was emerg- 
ing was the possible interpretation of nonextensive statis- 
tics as a kind of average of BG statistics, where the tem- 
perature (or scaled coupling constants) itself would be a 
stochastic variable. Such mathematical relations were in 
fact already present in the so called Hilhorst (see [64]), 
Prato [65] and Lenzi [66] formulae, however without any 
particular physical interpretation. 

The next step was recently done by Beck and Cohen 
when they proposed, through a Laplace transform, their 
superstatistics [67] . They generalized Eq. (84) into 



B{Ei) 



f 

Jo 



(86) 



where /(/3) is a quite general distribution ( d(}f{(3) = 
1). They consistently generalized Eq. (85) into 



1BC - 



(87) 



where (• • •) = Jq°° d(3{- ■ •)/(/?) (the subindex BC stands 
for Beck-Cohen and was introduced in [68] to avoid con 
fusion with q). Clearly, if /(/?) — fq{(3) we have that 

qsc = q- 

Of course, having a weight like that defined in Eq. 
(86) is necessary but not sufficient for having a statisti- 
cal mechanics. We also need an entropy functional, and 
a consistent manner of writing the constraints such as 
the energy and similar ones. The optimization of this 
entropy is expected to be extremal for the superstatisti- 
cal weight B{Ei). Such an entropy, from now on noted 
SgHPi}) {G stands for generalized), was in fact proposed 
in [68]. It has an imivocal relation with /(/?), i.e., one 
and only one functional Sq corresponds to each admissi- 
ble /(/?). Of course, /(/3) = /g(/3) imphcs and is implied 
by Sg = Sq. This sensibly general entropy Sq is concave 
by construction. It also happens to be experimentally 
robust, as proved in [69]. Such a convergence of two non- 
trivial properties (concavity and robustness) is certainly 
very satisfactory, and constitutes an argument favorable 
to considering Sq as a physical entropy (i.e., useful for 
thermal physics), and not just as one more theoretical 
information measure, among dozens that exist in the lit- 
erature of cybernetics, control theory, and other applied 
sciences. 

C. Connecting the differential-equation and the 
superstatistics paths 

The solution of Eq. (79) clearly corresponds to a par- 
ticular B{Ei), from now on noted Bqqi{Ei), and to its 
associated /(/3), from now on noted fqqi{(3). In partic- 
ular, we have that Bqi(Ei) equals the function given in 
Eq. (78) with {ai,aq,x) —^ (— — i?i). More partic- 
ularly. 



B2iiE) 



1 



l-f^(l-e/3i^) 



(88) 



If we normalize this function in order to get the probabil- 
ity distribution, take the limit 00, and multiply 
by the d = 3 photonic density of states, we recover, as 
already discussed, Planck's law. 

More details about Bqqi{Ei), fqq'{0), and their appli- 
cation to the recent high-precision experimental (by Bo- 
denschatz et al [70]) and computational (by Gotoh et al 
[71]) results for fully developed turbulence, can be found 
in [72]. Finally, the logical structure of the successive 
generalizations of classical BG statistical mechanics pre- 
sented in this paper is indicated in Fig. 2. The case of 
the Planck black-body distribution outside BG statisti- 
cal mechanics in the figure deserves a clarification. It is 
well known that, photons being bosons, this distribution 
naturally belongs to the quantum, version of BG statisti- 
cal mechanics, more precisely to Bose-Einstein statistics. 
To be more explicit, classical BG statistics is contained 
inside quantum BG statistics. Bose-Einstein statistics 
(as well as Fermi-Dirac statistics) is also contained in- 
side quantum BG statistics though outside from classi- 
cal BG statistics. All these distributions (classical BG, 
Bose-Einstein and Fermi-Dirac) can obviously be con- 
ceived as particular cases of the Beck-Cohen superstatis- 
tics since only a Laplace transform is involved (see Eq. 
(86)). What is special about Planck law is that it can 
be also conceived as a simple particular case of the spe- 
cial instance of superstatistics which simultaneously is a 
solution of Eq. (77). In other words, the celebrated dis- 
tribution can also be seen as a crossover between classical 
BG statistics (i.e., classical q — I statistics) and classical 
q = 2 statistics. Summarizing, it is in this specific sense 
that this distribution is, as indicated in Fig. 2 , outside 
classical BG statistics, and inside superstatistics. 



VI. CRITIQUES 



As the history of sciences plethorically shows to us, 
every possible substantial progress in the foundations 
of any science is accompanied by controversies. This 
is a common and convenient mechanism for new ideas 
to be checked and better understood by the scientific 
community. There is absolutely no reason to expect that 
statistical mechanics would be out of it. Quite on the 
contrary [73], given the undeniable fact that entropy is 
one among the most subtle and rich concepts in physics. 
Indeed, as c;loquc;ntly commented by Nicolis and Daems 
[73], "It is the strange privilege of statistical mechanics 
to stimulate and nourish passionate discussions related 
to its foundations [■■■]". We believe that some space 
dedicated here to such issues might well be useful at this 
stage. Therefore, let us address one by one some recent 
critiques that we are aware of. In fact they have all been 
replied in specific papers which we indicate here case 
by case for the interested reader. In what follows, for 
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economy of space, we shall restrict to a brief review of 
what we understand to be the main focus of each criti- 
cism, and what we believe to be the main reason for its 
rebuttal. 

(a) Vollmayer-Lee and Luijten critique: [74] 

Vollmayr-Lee and Luijten (VLL) presented in 2001 [74] 
a critique to nonextensivc statistical mechanics. They 
consider a Kac-potential approach of nonintegrable in- 
teractions. They consider a d-dimensional classical fluid 
with two-body interactions exhibiting a hard core as 
well as an attractive potential proportional to r^" with 
< a/d < 1 (logarithmic dependence for a/d = 1; 
VLL use the notation r = a). In their approach, they 
also include a Kac-like long-distance cutoff R such that 
no interactions exist for r > R, and then discuss the 
R^ 00 limit. They show that the exact solution within 
Boltzmann-Gibbs statistical mechanics is possible and 
that no surprise (see VLL Ref. [12] and references 
therein) - it exhibits a mean field criticality. Moreover, 
the authors argue that very similar considerations hold 
for lattice gases, 0{n) and Potts models. 

VLL state "Our findings imply that, contrary to some 
claims, Boltzmann-Gibbs statistics is sufficient for a 
standard description of this class of nonintegrable inter- 
actions. ", and also that "we show that nonintegrable in- 
teractions do not require the application of generalized 
q-statistics. ". In our opinion, these statements may mis- 
guide the reader. The critique was rebutted in [75], whose 
main points are summarized here. Indeed, the VLL dis- 
cussion, along traditional lines, of their specific Kac-like 
model only exhibits that Boltzmann-Gibbs statistical me- 
chanics is as more than one century of brilliant successes 
guarantees! - necessary for calculating, without doing 
time averages, a variety of thermal equilibrium proper- 
ties; by no means it proves that it is sufficient, as we shall 
soon clarify. Neither it proves that wider approaches such 
as, for instance, nonextensivc statistical mechanics (VLL 
Refs. [6,31] and present [12-14]), or any other similar 
formalism that might emerge, are not required or con- 
venient. The crucial point concerns time, a word that 
nowhere appears in the VLL paper. The key role of t 
has been emphasized in several occasions, for instance in 
Fig. 4 of VLL Ref. [31]. For integrable or short-range 
interactions (i.e., for a/d > 1), we expect that the t oo 
and N ^ oo limits are commutable in what concerns the 
equilibrium distribution p{E) , E being the total energy 
level associated with the macroscopic system. More pre- 
cisely, we expect naturally that (excepting for the density 
of states) 

p{E) = lim lim p{E; N; t) = lim lim p{E; N: t) 

oc exp[-£;/fcT] {r/d > 1) (89) 

if the system is in thermal equilibrium with a thermostat 
at temperature T. In contrast, the system is expected 
to behave in a more complex manner for nonintegrable 



(or long-range) interactions, i.e., for < a/d < 1. In 
this case, no generic reason seem to exist for the t oo 
and N ^ oo limits to be commutable, and consistently 
we expect not necessarily equal results. The simplest of 
these results (which is in fact the one to be associated 
with the VLL paper, although therein these two relevant 
limits and their ordering are not mentioned) is, as we 
shall soon further comment, 

lim lim p{E; N; t) oc exp[-{E/N)/{kT/N)] . (90) 

JV— >oo t— >(X) 

N = [N^-'^/'^ - a/d]/[l - a/d] has been introduced in 
order to stress the facts that generically 

(i) E is not extensive, i.e., is not proportional to N 
but is instead E oc NN [more precisely, E is extensive if 
a/d > 1 (see [76] and VLL Refs. [4,5]), and it is nonex- 
tensivc if < a/d < 1]; and 

(ii) T needs, in such calculation, to be rescaled (a fea- 
ture which is frequently absorbed in the literature by arti- 
ficially size-rescaling the coupling constants of the Hamil- 
tonian), in order to guarantee nontrivial finite equations 
of states. Of course, for a = 0, we have N = N, which 
recovers the traditional Mean Field scaling. 

B\it; depending on the initial conditions, which deter- 
mine the time evolution of the system if it is assumed 
isolated, quite different results can be obtained for the 
ordering Imit^oo^^T^N^oo p{E: N;t). This fact has been 
profusely detected and stressed in the related litera- 
ture (see, for instance, VLL Ref. [31], present Refs. 
[39,43-45,77-80] and references therein). Unfortimately, 
this important fact has been missed in the VLL critique. 
Such metastable states can by no means be described 
within BG statistical mechanics. Even more, as shown 
in Section IV, there is nowadays increasing evidence that 
they might be intimately related to nonextensivc statis- 
tical mechanics. In any case, it is plain that, for such 
long-range Hamiltonians, BG statistics is necessary but 
not suSicient, in contrast with the VLL statements. 

(b) Nauenberg critique: [81] 

Some line of critique concerns whether the zeroth prin- 
ciple of thermodynamics and thermometry are consistent 
with nonextensivc statistical mechanics. Such question- 
ing is by no means new: a couple of dozens of papers are 
available in the literature which address this important 
point. It has been recently raised once again, this time 
by Nauenberg [81]. He concludes, among many other 
critiques, that it is not possible to have thermalization 
between systems with different values of q. It appears 
to be exactly the opposite which is factually shown in 
[46]. His critique is rebutted in [46]. One of the cru- 
cial points that is unfortunately missed in [81], concerns 
discussion of "weak coupling" in Hamiltonian systems. 
Indeed; if we call c the coiipling constant associated with 
long range interactions (i.e., < a/d < 1), we have that 
limjv->(x) liiiic_>o cA'' = 0, whereas limc^o liinjv^oo c/V 
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diverges. No such anomaly exists for short-range inter- 
actions (i.e., Oijd > 1). Indeed, in this simpler case, we 

have that lim^v^oo linic^o cN = limc^o liniiv^oo ciV = 0. 
The nonuniform convergence that, for long-range inter- 
actions, exists at this level possibly is related to the 
concomitant nonuniform convergence associated with 
the t ^ oo and N ^ oo limits discussed previously in 
this paper. All these subtleties are not mentioned in [81]. 

(c) Luzzi, Vasconcellos and Ramos critique: [82] 

Another line of critique concerns the "physicality" of 
Sq (see [82]). Or whether it could exist a "physical" en- 
tropy different from Sbg- Since such issues appear to be 
of a rather discursive/philosophical nature, we prefer to 
put these critiques on slightly different, more objective, 
grounds. We prefer to ask, for instance, (i) whether Sq is 
useful in theoretical physics in a sense similar to that in 
which Sbg undoubtedly is useful; (ii) whether q necessar- 
ily is a fitting parameter, or whether it can be determined 
a priori, as it should if we wish the present theory to be 
a complete one; (iii) whether there is no other way of 
addressing the thermal physics of the anomalous systems 
addressed here, very specifically whether one could not 
do so by just using Sbg] (iv) whether Sq is special in 
some physical sense, or whether it is to be put on the 
same grounds as the twenty or thirty entropic functionals 
popular in cybernetics, control theory, and information 
theory generally speaking. Such questions have received 
answers in [83-88] and elsewhere, (i) The usefulness of 
this theory seems to be answered by the large amount of 
applications it has already received, and by the ubiquity 
of the (/-exponential form in nature, (ii) The a priori 
calculation of q from microscopic dynamics has been 
specifically illustrated in Section III. (iii) The optimiza- 
tion of Sq, as well as of almost any other entropic form, 
with a few constraints has been shown in [84] to be equiv- 
alent to the optimization of Sbg with an infinite number 
of appropriately chosen constraints. Therefore we could 
restrain to the use of Sbg if we absolutely wanted that, 
similarly to the fact that, instead of using the extremely 
convenient Keplerian ellipse for the planetary orbits, we 
could equivalently use an infinite number of Ptolemaic 
epicycles. It is however appreciably much simpler to rep- 
resent a complex structure of constraints into the single 
index q =/= 1 (in analogy with the fact that the ellipticity 
of a Keplerian orbit can be simply specified by a single 
parameter, namely the eccentricity of the ellipse), (iv) 
The entropy Sq shares with Sbg an impressive set of 
important properties, which namely includes concavity, 
stability and finiteness of the entropy production per unit 
time, Vg > 0. The difficulty of simultaneously satisfying 
all these three properties can be measured by the fact 
that the Renyi entropy (usefully used in the geometric 
characterization of multifractals) satisfies none of them 
for all g > 0. Such features point Sq as being very special 
for thermostatistical purposes. 



(d) Zanette and Montemurro critique I: [89] 

In a recent paper [89], Zanette and Montemurro re- 
analyze the molecular dynamics approach and results 
presented in [43] for the infinitely-long-range interacting 
planar rotators already discussed here. They especially 
focus on the time dependence of the temperature T{t) 
defined as the mean kinetic energy per particle. For to- 
tal energy slightly below the second-order critical point 
and a non-zero-measure class of initial conditions, a 
long-standing nonequilibrium state emerges before the 
system achieves the terminal BG thermal equilibrium. 
When T{t) is plotted, as done by virtually all authors, 
by using a linear scale for T and a logarithmic scale 
for t, an inflection point exists. If we call tcrossover 
the value of t at which the inflection point is located, 
it has been repeatedly verifled numerically by various 
authors, including Zanette and Montemurro [89], that 
lmijs[^octcrossover{l^) divcrgcs. Therefore, if the sys- 
tem is very large (in the limit N oo, mathematically 
speaking) it remains virtually for ever in the anomalous 
state, currently called by many authors, quasi- stationary 
state or metastable state. Zanette and Montemurro point 
out (correctly) that, if a linear scale is used for t, the 
inflection point disappears. From this, these authors 
conclude that this well known metastable state is but 
a kind of mathematical artifact, and no physically rele- 
vant quasi-stationarity exists. It is like if the high-to-low 
energies crossing occurring, at a given temperature, in 
Fermi-Dirac statistics had no physical meaning! Indeed, 
if instead of using the linear scale for the energies we 
were to use a faster scale (e.g., an exponential scale), the 
well known inflection point will surely disappear. Nev- 
ertheless, there is no point to conclude from this that 
the textbook crossing in Fermi-Dirac statistics is but a 
mathematical artifact. In fact, any inflection point on 
any curve will disappear by sufficiently "accelerating" 
the abscissa. Coming back to the system of rotators, 
what indeed appears to happen is that, for increasingly 
large N, T{t) remains constant, and different from the 
BG value, within an error bar which appears to vanish 
in the N oo limit. This effect appears in an even 
more pronounced way because of a slight minimum that 
T{t) presents just before going up to the BG value. This 
intriguing minimum had already been observed in [43] 
and has been detected with higher precision in [89]. A 
partial rebuttal of this critique is presented in [77] . 

(e) Zanette and Montemurro critique II: [90] 

Soon after the previous one, Zanette and Montemurro 
advanced a second critique [90] objecting the validity 
of nonextensive statistical mechanics for thermodynam- 
ical systems. This line of critique concerns the ubiquity 
of the (jT-cxponential form as a stable law in nature. 
The argument essentially goes that only Gaussians and 
Levy distributions would be admissible, because of the 
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respective central limit theorems. Such question has 
been addressed long ago in [55] (ignored in [90]), and 
once again more recently in [91], as a rebuttal to [90]. 
The answer basically reminds that the time stability of 
the distribution only has a mathematically well defined 
meaning if we also specify the time evolution or time 
composition law. The just mentioned central limit theo- 
rems only apply for convolution time evolution, certainly 
not for nonmarkovian evolutions as those illustrated in 
Section IV. 

(f) Zanette and Montemurro critique III: [92] 

Soon after the second critique, Zanette and Monte- 
murro advanced a third one [92]. This time the ob- 
jection addresses the rest of the systems, i.e., the non 
thermodynamical ones. It is argued by these authors 
that non thermodynamical applications of nonextensive 
statistics are ill-defined, essentially because of the fact 
that any probability distribution can be obtained from 
the nonextensive entropy Sg by conveniently adjusting 
the constraint used in the optimization. We argue here 
that, since it is well known to be so for any entropic 
form and, in particular, for the Boltzmann-Gibbs en- 
tropy Sbg, the critique brings absolutely no novelty to 
the area. In other words, it has nothing special to do 
with the nonextensive entropy Sq. In defense of the usual 
simple constraints, typically averages of the random vari- 
able Xi or of x| (where Xi is to be identified according 
to the nature of the system), we argue, and this for all 
entropic forms, that they can hardly be considered as 
arbitrary, as Zanette and Montemurro seem to consider. 
Indeed, once the natural variables of the system have 
been identified, the variable itself and its square obvi- 
ously are the most basic quantities to be constrained. 
Such constraints are used in hundreds (perhaps thou- 
sands) of useful applications outside (and also inside) 
thermodynamical systems, along the information theory 
lines of Jaynes and Shannon, and more recently of A. 
Plastino and others. And this is so for Sbg, Sq, and any 
other entropic form. If, however, other quantities are 
constrained (e.g., an average of .t'^ or of [.t]"^) for specific 
applications, it is clear that, at the present state-of-the- 
art of information theory, and for all entropic forms, this 
must be discussed case by case. The full rebuttal of this 
critique can be found in [93]. 

(g) About ordinary differential equations: 

The remarks in Sections LA and II. A related to ordi- 
nary differential equations might surprise some readers, 
hence deserve some clarification. Indeed, in virtually all 

the textbooks of statistical mechanics, fimctions such as 
the energy distribution at thermal equilibrium are dis- 
cussed using a variational principle, namely referring to 
the entropy functional, and not using ordinary differen- 
tial equations and their solutions. In our opinion, it is 
so not because of some basic (and unknown) principle of 



exclusivity, but rather because the first-principle dynam- 
ical origin of the BG factor still remains, mathematically 
speaking, at the status of a dogma [94]. Indeed, to the 
best of our knowledge, no theorem yet exists which es- 
tablishes the necessary and sufficient conditions for being 
valid the use of the celebrated BG factor. Nevertheless, 
one must not forget that it was precisely through a dif- 
ferential equation that Planck heuristically found, as de- 
scribed in his famous October 1900 paper [2], the black- 
body radiation law. It was only in his equally famous 
December 1900 paper that he made the junction with 
the Boltzmann factor by assuming the — at the time, 
totally bizarre — hypothesis of discretized energies. 

A further point which deserves clarification is why have 
we also interpreted the linear ordinary differential equa- 
tion in Section LA as providing the typical time evolution 
of both the sensitivity to the initial conditions and the 
relaxation of relevant quantities. Although the bridging 
was initiated by Krylov [1] , the situation still is far from 
completely clear on mathematical grounds. However, 
intuitively speaking, it seems quite natural to think that 
the sensitivity to the initial conditions is precisely what 
makes the system to relax to equilibrium, and therefore 
opens the door for the BG factor to be valid. In any case, 
although some of the statements in Section I. A are (yet) 
not proved, this by no means implies that they are gener- 
ically false. Furthermore, they provide what we believe 
to be a powerful metaphor for generalizing the whole 
scheme into the nonlinear ordinary differential equations 
discussed in Section II. A. Interestingly enough, the q- 
exponential functions thus obtained have indeed proved 
to be the correct answers for a sensible variety of specific 
situations reviewed in the rest of the present paper, and 
this for all three interpretations as energy distribution 
for the stationary state, time evolution of the sensitivity 
to the initial conditions, and time evolution of basic re- 
laxation functions. 



VII. FINAL COMMENTS 

We all know that the concepts of energy and entropy 
are cornerstones of contemporary physics. They are both 
at the heart of thermodynamics, a set of connections 
that regulate the laws of the macroscopic world. Sta- 
tistical mechanics bridges microscopic dynamics (Newto- 
nian, quantum, relativistic, or any other) with the macro- 
scopic behaviors we observe in nature. Boltzmann-Gibbs 
statistical mechanics constitutes an impressive illustra- 
tion of how the logics of these micro-macro connections 
is structured. This theory is based on the BG expres- 
sion Sbg for the entropy. Following along the lines of 
Einstein 1910 and of many others, it is our belief that 
the specific mathematical form that Sbg has must ulti- 
mately descend from microscopic dynamics. More pre- 
cisely, this connection is known to rely on quick mixing 
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in phase space, ergodicity, and related nonlinear dynam- 
ical concepts. This is the statistical mechanics of short 
range interactions, short-range microscopic memory, ul- 
timately leading to a simple, uniform occupancy of the 
allowed phase space (the "equal probability" hypothesis 
for the microcanonical ensemble). The question arises 
naturally: what happens when the interactions are long- 
ranged, when the microscopic memory is heavily non- 
markovian, when the geometry of occupation of phase 
space is complex, (multi) fractal for instance, or some 
other hierarchical structure? Can we still devise theoret- 
ical techniques similar to those of BG statistical mechan- 
ics, which would once again bridge with the laws of the 
macroscopic world? We believe that we can. A central 
point therefore is: what mathematical expression should 
or could we use for the entropy? How can we adequately 
generalize Sbg? 

The answer to such question obviously is far from triv- 
ial. Indeed, even in the framework of BG statistical me- 
chanics, we do not know how to rigorously make the 
expression of Sbg descend from microscopic dynamics, 
plus possibly some other generic logical requests. In some 
sense, the best we know nowadays about such fundamen- 
tal question lies at the level of the necessary and sufficient 
conditions proposed by Shannon, by Khinchin, and sim- 
ilar constructs. But they do not start from dynamics. 

Since the direct connection between normal (ergodic) 
microscopic dynamics and Sbg is yet not totally clear, 
there is no surprise that the same happens for anomalous 
microscopic dynamics. This is the fundamental reason 
for which we have proceeded to generalize Sbg into Sq 
through a metaphoric path, as illustrated in Section II. 
The (/-generalizations of the Shannon and the Khinchin 
theorems have already been established by Santos and 
by Abe respectively. A variety of arguments qualify the 
conditions under which Sg is unique (see [83] and refer- 
ences therein). Still, how to make Sg (or even Sbg^ ^is 
we said before) descend from microscopic dynamics re- 
mains an open question. Consistently remains an open 
question what exactly have in common all the physical 
cases indicated in Section III, and what may exactly be 
the geometrical striicture which reflects the dynamical 
occupancy of the phase space in such cases. We have 
conjectured that it may well be similar to a scale- free ge- 
ometry like that of the Albert-Barabasi model, but this 
remains to be proved. 

Another intriguing question is the connection with 
aging, a property that long-range-interacting Hamilto- 
nians share with many glasses, spin-glasses, and other 
metastable systems intensively studied nowadays. Such 
long-range Hamiltonians satisfy, even in the longstand- 
ing aging phase (whose duration diverges with N), a 
zeroth principle of thermodynamics, which concretely 
opens the path to thermodynamics. Furthermore, a sen- 
sible amount of connections with nonextensive statistical 
mechanics are already available in the literature. How- 
ever, the exact dependence of q on {a,d), possibly on 
a/d, is still to be unambiguously established. Many re- 



searchers around the world are presently working on that 
fascinating problem. 

Finally, we should emphasize that everything that 
we know today neatly points towards the scenario that 
nonextensive statistical mechanics satisfies the 0-th, first, 
second and third principle of thermodynamics (see, for 
instance, [95]). This is deeply interesting, since that im- 
plies that the basic laws of thermodynamics are stronger 
than the role reserved for them within Boltzmann-Gibbs 
statistical mechanics. Further analysis of the foundations 
and thermodynamical connections of statistical mechan- 
ics in general, and of nonextensive statistical mechanics 
in particular, are certainly fascinating and very welcome. 
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FIG. 1. Successive versions for the canonical-ensemble en- 
ergy-constraint within nonextensive statistical mechanics re- 
garding properties (i-vii) indicated in Section III assuming 
fixed the Lagrange parameter 13. In the 1988 paper [12] two 
possible constraints for the energy were indicated , namely 
^iPiEi = constant and '^iP'iEi — constant, but only the 
former was developed therein. In the 1991 paper [13], the 
second form for the constraint was developed, which enabled 
the connection with thermodynamics. In the 1998 paper [14], 
the constraint using escort distributions was adopted and de- 
veloped, i.e., Xli Pi -^«/ Xii Pi ~ constant. This last form pro- 
vides for the points (i-vii) the same answers as within BG 
statistical mechanics, and is presently believed to be its cor- 
rect g-generalization. 



FIG. 2. Logical structure of the successive generaliza- 
tions of classical BG statistical mechanics presented in 
this paper, with indication of their paradigmatic differen- 
tial equations, and their /(/3) superstatistical distributions; 
/i(/3) = /sg(/S) ;,f<j(/3) = -distribution. The correspond- 
ing stationary state probability distributions and entropies 
are given by Eqs. (9) and (10), for BG statistical mechanics, 
by Eqs. (30) and (34), for nonextensive statistical mechan- 
ics, and by Eq. (86) and So given in [68,69], for generalized 
nonextensive statistical mechanics. The q ^ q crossover sta- 
tistical mechanics includes as a particular case the g 1 
crossover one, which corresponds to Eq. (88). If we take 
into account in this equation the normalization factor, the 
132/ Pi — > 00 limit, and the photonic density of states, we re- 
cover Planck's law for the black-body radiation. See the text 
in Section V.C for clarification of this case. 



22 



